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Abstract 

It is shown that geodesic motion on the GL{n, M) group manifold endowed with 
the bi-invariant metric ds'^ = ti (^g~^dg^'^ corresponds to a generahzation of the 
hyperbohc n-particle Calogero-Moser-Sutherland model. In particular, considering 
the motion on Principal orbit stratum of the SO{n,M) group action, we arrive at 
dynamics of a generalized n-particle Calogero-Moser-Sutherland system with two 
types of internal degrees of freedom obeying SO{n,U.) ^ SO{n,M.) algebra. For the 
Singular orbit strata of SO{n,M) group action the geodesic motion corresponds 
to certain deformations of the Calogero-Moser-Sutherland model in a sense of de- 
scription of particles with different masses. The mass ratios depend on the type 
of Singular orbit stratum and are determined by its degeneracy. Using reduction 
due to discrete and continuous symmetries of the system a relation to Il^n Euler- 
Calogero-Moser-Sutherland model is demonstrated. 
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1 Introduction 



Almost twenty years ago a possibility was discovered to maintain the integrability of 
Calogero-Moser-Sutherland models [1] (classification and description can be found in [2]) 
supposing that the particles moving on a line have additional internal degrees of freedom 
[3,4]. Latter it was shown [5,6] that the generic elliptic Calogero-Moser-Sutherland type 
system, which consists of n-particles on a line interacting with pairwise potential in the 
form of Weierstrass elliptic function V{z) = p{z), admits the following generalization 

n 1 

H = ^Y.P^ + 2^ - Xj) . (1) 

1=1 i^j 

Here apart from the canonical pairs {xi,pi), describing the position and the momenta of 
particles and obeying nonvanishing Poisson brackets 

{xi,pj} = 6ij, (2) 

the "internal" degrees of freedom fab which satisfy the algebra 

{fab, fed} = ^bcfad — ^adfcb (3) 

are included. The recent comprehensive discussion of the integrability of generic spin 
Calogero-Moser-Sutherland systems can be found in the papers [7,8]. One of the most 
effective and transparent way to convince in the integrability of a Hamiltonian system 
is to find a known higher- dimensional exactly solvable model, whose dynamics on a cer- 
tain invariant submanifold coincides with the dynamics of the given Hamiltonian system. 
This method is known as symplectic reduction method [9-11]. The Calogero-Sutherland- 
Moser systems with the so-called degenerate cases of the potential, when the Weierstrass 
elliptic function p{z) reduces to l/sinh^2;, 1/sin^z or to the rational function 
have a well-known interpretion as symplectic reductions of geodesic motions on sym- 
metric spaces [12,13]. Furthermore, it was argued in [14] that symplectic reduction relates 
the elliptic Calogero-Moser-Sutherland system with certain integrable Hamiltonian sys- 
tem on the cotangent bundle to the central extension of two-dimensional Lie algebra of 
SL{n, C)-valued currents on the some elliptic curve. New types of generalizations of the 
spin Calogero-Moser-Sutherland systems with nonstandard spin interactions have been 
constructed in [15] using discrete symmetries of the model. 

In the present Letter we shall exploit the idea of symplectic reduction considering certain 
generalization of the Calogero-Sutherland-Moser model. Namely, we shall consider the 
integrable finite-dimensional model corresponding to the geodesic motion on the general 
linear matrix group with a positive determinant GL^{n, M),[^ endowed with the left- and 
right-invariant metric ds"^ = tr (g~^dg)'^, where g G GL(n, R). In terms of this bi-invariant 

^ Hereafter we shall omit the upper index + to simplify the expressions. 
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metric on GL(n,M.) group manifold the equations of motion for the corresponding dy- 
namical system are encoded in the Lagrangian [9,10] 



LGL = ltT{g''gy , (4) 

where overdot denotes differentiation with respect to time. Bellow we shall represent 
the Hamiltonian corresponding to Lagrangian (4) in terms of a special parameterization, 
adapted to the action of SO{n,M.) symmetry group of the system. We shall demonstrate 
that on the Principal orbit stratum of S'0(n, M) group action the resulting Hamiltonian 
defines a new generalization of the Calogero-Sutherland-Moser model by introducing two 
internal variables "spin" and "isospin" . Furthermore, performing the Hamiltonian reduc- 
tion owing to two types of symmetry: continuous and discrete, we show how to arrive at 
the conventional Hamiltonian of Euler-Calogero-Sutherland model 

-in 1 " 

H-IEpI + IE-^ ^ (5) 

2 ^ 8 ^ smh [Xi - Xj) 

with internal variables lab — — ha, obeying the 5*0(71, M) Poisson bracket algebra 

{lab, led} — ^achd — ^adhc + hdlac — ^beLd ■ (6) 



Another interesting systems arise when the dynamics takes place on the Singular orbit 
strata of the SO{n,M.) group action. We found in this case new models representing a 
certain class of mass deformed Calogero-Moser-Sutherland models. In particular, for the 
case of GL{3,M) group, our analysis shows that the dynamics on Singular orbit stratum 
with isotropy group 50(2) Z2, corresponds to Calogero-Moser-Sutherland model, de- 
scribing two particles, whose mass ratios is mi : m2 = 2 : 1 (see eq. (68)). The question of 
integrabihty of the mass deformed Calogero-Moser-Sutherland models has been discussed 
in [16] and references therein. 



2 Geodesic motion on the Principal orbit stratum 



2.1 Symmetries and dynamics 



If we choose the elements of the matrix g e GL{n,M.) as Lagrangian coordinates, the 
Euler- Lagrange equations obtained from the Lagrangian (4) can be represented in the 
form of current conservation 

^M)-0. (7) 
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This form allows to find the general solution, i.e. the geodesies of the bi-invariant metric 
are given by 

g{t)=giO)exp{tJ), (8) 

where g{0) and J arc two arbitrary constant matrices. The special choice of these matri- 
ces corresponds to the particular solutions describing the motion on a certain invariant 
submanifold. 



To show that (7) are equations of motion of many-particle system representing a certain 
generalization of the Calogero-Moser-Sutherland model, it is useful to pass to the Hamilto- 
nian form of the geodesic motion on GL{n, M) group. Performing Legendre transformation 
of the Lagrangian (4) 

T 9 Lgl ( -1 



we arrive at the canonical Hamiltonian 



2 

generating the Hamilton equations of motion 



i^GL = ^tr(7r^5)' (10) 



9 = {9,Hgl} = (11) 

TT = {tT ,Hgl} = -T!- g^ TT . (12) 

The nonvanishing Poisson brackets between the fundamental phase space variables {gab, T^ab) 
are 

{gab^T^cd} ^ SacSbd- (13) 



Prom now on the purpose of the present paper will be to rewrite this Hamiltonian in terms 
of coordinates, adapted to the symmetry possessing the system. At first we would like to 
analyze the following symmetry action of the SO{n, R) group on GL{n, R) 

g^g' = Rg (14) 

with time-independent orthogonal matrix R. In order to consider the configuration space 
as manifold with orbit and slice structure with respect to this action, it is convenient to 
use the polar decomposition [17] for an arbitrary element of the GL{n, M) group. For the 
sake of technical simplicity we investigate in details the G'L(3,R) group hereinafter, i.e. 

9 -OS, (15) 
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where S is a, positive definite 3x3 symmetric matrix, and O(0i, 02, <t>3) = e'^^'^'^e'^^'^^e'^^'^^ 
is an orthogonal matrix with 5*0(3, M) generators in adjoint representation {Ja)ij = £iaj- 
Since the matrix g represents an element of GL(3, K.) group, we can treat the polar de- 
composition (15) as a uniquely invertible transformation from the configuration variables 
g to a new set of Lagrangian variables: six coordinates Sij and three coordinates (pi. In 
terms of these new variables the Lagrangian (4) can be rewritten as 

LoL^ltr (Ql + SS-'Y , (16) 

where 0^ := O is a left-invariant 1-form on the 5*0(3, M) group. To find the cor- 
responding Hamiltonian we note that the polar decomposition (15) induces the point 
canonical transformation from variables {gab, T^ab) to new canonical pairs {Sab, Pab) a-nd 
(0a, Pa) obeying the nonvanisning Poisson bracket relations 

{Sab , Pcd} = 2 i^ac hd + ^ad ^bc) , (17) 
{(t>a,Pb] = 5ab- (18) 

The expression of the old TTab as a function of the new coordinates is 

TT = O (P - kaJa) , (19) 

where 

K = 1^^ (Vb - ^bmn {SP)^^) , (20) 

Jik — Sik — SiktrS and rja are three left-invariant vector fields on 5'0(3,M) group 

Vi = Pi + cos 03 P2 - cot 02 sin 03 P3 , (21) 

sm 02 

r COS 03 

ri2 = ■ , -Pi - sin 03 P2 - cot 02 cos 03 P3 , (22) 
sm 02 

V3-P3. (23) 



Hence, in terms of the new variables, the canonical Hamiltonian (10) takes the form 

Hgl = ^ tr {PSf + ^ tr {JaSJbS) kah , (24) 

where the canonical variables {Sab, Pab) are invariant under the transformation (14), while 
the angular variables (0a, Pa) undergo changes generating by the right-invariant KiUing 
vector fields 77^ 
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rji^- sin 01 cot 02 Pi + cos 0i P2 + , (25) 

77^= COS 01 cot 02 Pi + sin 01 P2 - ^^^^^3, (26) 

sm 02 

Vs = ^'i , (27) 

whose Poisson brackets with the left-invariant vector fields r)^ vanish {rj^ , rj^} — 0. 

Now we pass to analysis of another type of symmetry action of the orthogonal group. The 
Lagrangian (4) is invariant under the transformations 

g^g'^R^gR (28) 

with constant orthogonal matrix R G 5*0(3, M). After implementation of the polar de- 
composition the symmetry transformation reads 

S'^R^SR, O'^R^OR. (29) 

The orbit space of the action S ^ R^ S R of the 5*0(3, M) group in the space of 3 x 3 
symmetric matrices S is given as a quotient 5/5*0(3, M). The quotient space 5/5*0(3, M) 
is a stratified manifold; orbits with the same isotropy group are collected into strata and 
uniquely parameterized by the set of ordered eigenvalues of the matrix S : xi < X2 < x^. 
The strata are classified according to the isotropy groups which are determined by the 
degeneracies of the matrix eigenvalues: 

(1) Principal orbit-type stratum, when all eigenvalues are unequal Xi < X2 < Xs, with 
the smallest isotropy group Z2 C?> ^2 • 

(2) Singular orbit-type strata forming the boundaries of the orbit space with 

(a) two coinciding eigenvalues (e.g. Xi = X2), when the isotropy group is 5*0(2) 0^2 . 

(b) all three eigenvalues are equal {xi — X2 = x^), here the isotropy group coincides 
with the isometry group 5'0(3, M). 

To write down the Hamiltonian describing the motion on the Principal orbit stratum, we 
introduce coordinates along the slices x and along the orbits Namely, since the matrix 
S is positive definite and symmetric, we use the main-axes decomposition in the form 

S = R^ix) e'"" R{x) , (30) 

where R{x) G 5*0(3, M) is an orthogonal matrix parameterized by three Euler angles 
X = {Xi: X2: Xs), and the matrix e^^ is diagonal e^^ = diag ||e^^i, e^^^ g2a;3|| ^j^g momenta 
Pi and p^., canonically conjugated to the eigenvalues Xi and the angles Xi correspondingly, 

{xi ,pj}^ Sij , {xi , PxJ = ^ij > (31) 
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can be found using the condition of canonical invariance of the symplectic 1-form 

3 3 3 

J2 Pij Sij dt = J2Pi^idt + Yl Pxi Xi dt . (32) 

i,j=l j=l i=l 

The original momenta are expressed in terms of the new canonical pairs {xi,pi) and 
iXhPxi) as 



\o=l a=l J 



(33) 



with 



Pa ^ I Pa: (34) 



2 



Va — TT^ 7 ) (cychc permutation a^h^ c) . (35) 

4 sinh(x6 — Xc) 

In the representation (33), we introduced the orthogonal basis for the symmetric 3x3 
matrices a a — (eta, eta), a = 1, 2, 3 with the scalar product 

tr(a„ Qib) = , tr(Q;a «(,) = 25ab , tT{aa ab) = (36) 

and the 5'0(3, M) right-invariant Killing vectors 



^1 = - sm xi cot X2 Pxi + cos Xi Px2 + Pxz > (37) 

smx2 

COS /^i 

^2 = COS xi cot X2 Pxi + sin Xi Px2 - Px3 > (38) 

sinX2 

^3 = Pxi ■ (39) 

Thus, after passing to main-axes variables {xi,Pi) and (Xi)Pxi)' canonical Hamiltonian 
reads 

H^^-^± f I ^ ' E ^ ^ • (40) 

8„^i 16^^^^)Sinh2(x6-Xc) 4(^^^^ cosh2(x6-Xc) 

Here (a6c) means cyclic permutations a ^ b ^ c. Hence we conclude that the integrable 
dynamical system describing a free motion on the Principal orbit stratum can be inter- 
preted in the adapted basis, as Generalized Euler-Calogero-Moser-Sutherland model. The 
generalization consists in the introduction of two types of internal dynamical variables ^ 
and 7] — "spin" and "isospin" degrees of freedom. From their explicit expressions (see eqs. 
(21)-(23) and (37)-(39)) it follows that they satisfy 50(3, M) 50(3, M) Poisson bracket 
algebra 
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{Va,Vb} = -eabcVc, (41) 
, en = Babc^ , (42) 
{^a , } = . (43) 

2.2 Lax-pair for the Generalized Euler-Calogero-Moser-Sutherland model 

In order to find a Lax representation for tlie generalized Euler-Calogero-Moser-Sutlierland 
model (40) let us consider the integrals of the geodesic motion on the Principal orbit 
sratum. The integrals of motion can be written in Hamiltonian form, following from (7), 
as 

Jab = {T^'^g)ab ■ (44) 

The algebra of this integrals realizes on the symplectic level the GL(n, R) algebra 

{Jab-, Jed} = hcJad — SadJcb ■ (45) 

After the transformation to scalar and rotational variables (30), the expression for the 
current J reads 

1 ^ 

a=l 



where 



and 



= E ^ ef coth(xfc - Xc) + (^RamV^ + ^ ) tanh(a;b - Xc) (47) 

(abc) 



Ja = RamV^ + ^a- (4^ 



Using expressions (46) for the integrals Jab the classical equations of motion for General- 
ized Euler-Calogero-Moser-Sutherland model can be rewritten in the Lax form[^ 

L=[A,L], (49) 

where the 3x3 matrices are given explicitly as 



We set iiere aside the constructions of the Lax pairs with a spectral parameter. The Lax repre- 
sentations with a spectral parameter for the spin Calogero-Moser-Sutherland models associated 
with the root systems of simple Lie algebras were constructed in [7,8,18]. 
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^ Pi L+^L^"^ 



^ -A3, A2 ^ 

A3, 0, -Ai 
\^-A2, Ai, J 



(50) 



Entries and of the matrices (50) are given as 



and 



_ So ^'-am'lm ^ 2 

" 2 sinh^(a;6 — Xc) cosh^(a;fe — Xc) ' 
where (a, b, c) means cychc permutations of (1, 2, 3). 



(52) 



Below relations to the standard Euler-Calogero-Moser-Sutherland model (5) will be demon- 
strated. 



2.3 Reduction to Euler-Calogero-Moser-Sutherland model 



2.3.1 Reduction using discrete symmetries 

Now we shall demonstrate how the IIA3 Euler-Calogero-Moser-Sutherland model arises 
from the canonical Hamiltonian (10) after projection onto a certain invariant submani- 
fold determined by discrete symmetries. Let us impose the condition of symmetry of the 
matrices g e GL(3,M) 

^^'^=^-/ = 0. (53) 

In order to find an invariant submanifold, it is necessary to supplement the constraints 
(53) with the new ones 

^(2) ^ - TT^ = . (54) 

Indeed, using the Hamilton equations (11) and (12), one can check that the surface defined 
by the set of constraints a = {ip^^\ ip^"^^) represents an invariant submanifold in the 
GL(3,M) phase space 

^^^'^1*^=0 = 0, V'^')|*^=o = (55) 
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and the dynamics of the corresponding induced system is governed by the reduced Hamil- 
tonian 

1 2 

^GL(3,R)kA=0 = 2 ■ (56) 

The matrices g and tt are now symmetric nondegenerate matrices, and one can be con- 
vinced that this expression leads to the Hamiltonian of the IM3 Euler-Calogero-Moser- 
Sutherland model. To prove this, it is necessary to note that the Poisson matrix Cab = 
\\{ip^^\'4'^'^^}\\ is not degenerate and after projection on the invariant submanifold, the 
canonical Poisson structure is changed according to the Dirac prescription 

{F, G}n = {F, G}pB - {F, MC^Mb, G} (57) 

for arbitrary functions F and G. The resulting fundamental Dirac brackets are 

{gab , T^cdjo = ^ {Sac hd + ^ad he) ■ (58) 

If now we introduce as in Section (2.1) the main-axes variables (xq, Pa) and {xa, Pxa) 
instead of the symmetric variables {gah)T^ab) and use the representations (30) and (33), 
then one can convince that the projected Hamiltonian (56), up to the time rescaling 
t ^ 4t, governs the same dynamics as the Hamiltonian (5) of \1A^ Euler-Calogero-Moser- 
Sutherland model with the intrinsic spin variables lij — SijkCk- 



2.3.2 Reduction due to continuous symmetry 

Let us now derive a reduced Hamiltonian system employing certain continuous symmetries 
of the model. For the Hamiltonian (10) all angular variables are gathered in the three left- 
invariant vector fields rj^ and thus the corresponding right-invariant fields rj^ (25)-(27) 
are integrals of motion 

{77„^, Hgl} = . (59) 
The surface in the phase space determined by the constraints 

= (60) 

defines an invariant submanifold. These constraints obey the algebra {q^^ rj^} = EabcV^ 
and according to the Dirac terminology [19,20] are first class constraints, this means 
that after projection on the constraint shell (60), the corresponding cyclic coordinates 
disappear from the projected Hamiltonian. To prove this one can use the relation r)^ — 
OabVb between the left and the right-invariant Killing vector fields. Then, after projection 



10 



to the constraint surface (60), the Hamiltonian (40) reduces to 



1 ^ 1 

8V 4^^^smh 2{xb-Xc) 

After rescahng of the variables 2xa i— > Xa, one can be convinced that the derived Hamil- 
tonian coincides with the Euler-Calogero-Moser-Sutherland Hamiltonian (5), where the 
intrinsic spin variables are — SijkCk- 

As it was outlined in Section (2.2) apart from the integrals rj^ the system (10) possesses 
the integrals (46). Using these integrals one can choose different invariant submanifold 
and to derive the corresponding reduced system. Here we would like only to mention that 
after performing reduction to the surface defined by the vanishing integrals ja = 0, we 
again arrive at the Euler-Calogero-Moser-Sutherland system. 



3 Geodesic motion on the Singular orbit strata 



In the previous sections we have investigated the geodesic motion on the Principal orbit 
stratum, i.e. under the supposition that the symmetric matrix S in the polar representa- 
tion (15) has three different eigenvalues. We now turn our attention to dynamical system 
corresponding to the geodesic motion on the Singular orbit strata. For the sake of techni- 
cal simplicity, we restrict ourselves to invariant submanifold of the phase space, defined by 
7]^ = and consider a geodesic motion on the Singular orbit stratum with two coinciding 
eigenvalues of the matrix S in the case of GL(3,M) group. Below we use two alternative 
methods. At first a special parameterizations of the subspace of 3 x 3 symmetric matrices 
with two coinciding eigenvalues are exploited and as a result we find that the Hamiltonian 
system describing the geodesic motion on the 4-dimensional Singular orbit stratum is IIA2 
Calogero-Moser-Sutherland model with particle mass ratio nii : 1712 = 2:1. Afterwards, 
based on the observation that the Singular orbits of the configuration space represent 
the boundary of Principle orbit, using appropriate limiting procedure we derive from the 
Hamiltonian (40) again a mass deformed IIA2 Calogero-Moser-Sutherland model with the 
same particle mass ratio mi : m2 = 2:1. 



3. 1 Mass deformed Calogero-Moser-Sutherland model via explicit parameterizations of 
the Singular orbit stratum 



The Singular orbits have continuous isotropy groups and this leads to the modification 
of geodesic motion. For the case we are inteesting in, GL(3, M) group and two equal 
eigenvalues of the symmetric matrix S, it is SO (2) (g) Z2. 
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The linear space of the real symmetric nxn matrices with two coinciding eigenvalues has 
a real dimension [21] 

dim S{n) - dim S{2) + 1. (62) 

Hence, we are able to parameterize such a subspace of GL(3, M) group by four real inde- 
pendent parameters 

Sab = 5ab - 2 e"+^ sinh(x - y) UaU^ , (63) 
where Ua is a unit 3-dimensional vector 

ria = {sin 9 sin0, sin 6* cos0, cos 6*) . (64) 



We infer from the expression for the bi-invariant metric on the GL{n,M.) group that the 
metric induced on the 4-dimensional Singular orbit stratum parameterizing according to 
(63) is 

tr (^-^ dsf = 8 dx^ + 4dy'^ + 8 smh'^{x - y) [dO^ + sin 6 d^)"^) . (65) 

Therefore the Lagrangian L = | tr {s~^S^ on the Singular orbit stratum can be written 
as 

L = 4i;2 + 2y2 + 4sinh^(x-7/)n^ (66) 

where 

n' = ^' + sin2^(/.2. (67) 
The Legendre transformation gives the canonical Hamiltonian 

HoLim = r6P^ + -sPy+ I6sinh^(.-,) ' ^''^ 

where 

Now taking into account that is a constant of motion we convince that the geodesic 
motion with respect to the bi-invariant metric on the 4-dimensional Singular orbit stratum 
corresponds to 2-particle mass-deformed Calogero-Moser-Sutherland model with particle 
mass ratio mi : m2 = 2:1. Following this interpretation of the Hamiltonian system (68) in 
terms of particles, one can say that the motion on this Singular orbit stratum corresponds 
to some "gluing" of two particles and formation of a bound particle with double mass. It is 
apparent that using the center-mass coordinates the obtained 4-dimensional Hamiltonian 
system (68) can be reduced to 2-dimensional integrable model. 
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3.2 Mass deformed Calogero-Moser- Sutherland model via limiting procedure from the 
free motion on the Principal orbit stratum 



As it was mentioned above the Singular orbits with two coinciding eigenvalues form the 
boundary of the Principle orbit stratum. Based on this observation, now we would like 
to extend the Hamiltonian (40), given on the Principal orbit startum, to its boundary 
by introducing the constraints that force the dynamics specially to the neighborhood 
of the boundary and then use a limiting procedure. Since the Hamiltonian (40) has a 
singularities along the Singular orbits, we impose the following constraints on the phase 
space variables 

with positive small parameter e ^ 1 and constant a > 1. When e goes to zero the system 
tends to the Singular orbit stratum with two coinciding eigenvalues. Let us at first find 
the dynamical consistence of the constraints (70) considering instead of the Hamiltonian 
(40) the modified Hamiltonian Q 

H' = lj:pl + ^Y.^aVix,-x,)+u{x,-X2-e) + X (e3-e^+") , (71) 

a=l (abc) 

with Lagrangian multipliers u and A and V{x) = sinh~^(a;). The conservation in time of 
the constraints x^^"* and (p^^^ leads to the new constraints 

X^'^=Pi-P2, 0(^^=6, (72) 
while from the maintenance of (72), the Lagrangian multipliers can be found 



u 



16 
1 



^ ^1 V'ix2 -X3)-l ^2 V'{X^ - Xs) - V'{X^ - X2] 



A = - ^3 [V{X2 - X3) - V{Xi - X2)] . 



(73) 
(74) 



Hence we conclude that the constraint surface defined by (70) and (72) represents an 
invariant submanifold for the Hamiltonian (71). Because these constraints are second 
class in the Dirac terminology [11,19,20] we are able to replace the Poisson brackets by 
the Dirac ones according to (57) and let the constraint functions to vanish. One can easy 
verify that for the canonical variables {x,p) the corresponding nonzero fundamental Dirac 
brackets are 

{xi,pj}D = ^, i,j = 1,2, {x3,P3}d = 1, (75) 
^ Here we again restrict consideration by the case rj^ = 
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while for the angular variables we have 



{^a:^b}D = 0, a, 6 =1,2, 3. (76) 

Projecting the Hamiltonian H' to the constraint shell and then taking the limit e — > we 
obtain 

:= = \pI + IpI + r—r^ r ■ (77) 

'^0 4 8 16 smh (xi — Xsj 



The quantity in (77) is a constant of motion, that is the reminiscent of conserved total 
momentum + + ^| for the Euler-Calogero-Moser-Sutherland Hamiltonian (61). So, 
using the appropriate limiting procedure from the Principle orbit stratum we arrive at 
mass deformed Calogero-Moser-Sutherland model corresponding to the Singular stratum, 
labeled by two coinciding eigenvalues of 3 x 3 symmetric matrix. Finally, we establish 
a relation between (77) and the Hamiltonian (68) derived before. To achieve this it is 
necessary to rescale the momentum pi i— > 2~^pi so that variables (3:^1,^1,2^3,^3) obey 
canonical Poisson bracket relations instead of the Dirac brackets (75), identify variables 
X := Xi, y := 2:3, and constants P = ^f. As a result we arrive at Hamiltonian system 
which coincides with the mass deformed IM3 Euler-Calogero-Moser-Sutherland model, 
derived in the previous section using explicit parameterizations of the induced metric on 
Singular orbit stratum. 



4 Concluding Renicirks 



Nowadays we have revival of the interest to matrix models (see e.g. [22]) connected with 
the search of relations between the supersymmetric Yang-Mills theory and integrable 
systems (for a modern review see [23]). As it has been shown in the recent paper [24] 
the Euler-Calogero-Moser-Sutherland model with certain external potential describes the 
gauge invariant long- wavelength approximation of the SU{2) Yang-Mills field theory [25]. 
In the context of the consideration of higher dimensional gauge groups it is interesting to 
explore the mechanics on the general linear group manifold. In the present Letter we have 
considered the simplest version of geodesic motion on GL{n,'R) group manifold and and 
analyze the dynamics in the context of isometrics of the bi-invariant metric. Namely we 
intensively exploit the slice structure of GL{n, M) based on the existence of the Principal 
orbit stratum and Singular orbit strata of the SO{n,M) group action. We demonstrated 
that the free motion on the Principal orbit stratum corresponds to the integrable many- 
body system of free particles on a line with two types of internal variables called "spin" 
and "isospin", which is a generalization of the Euler-Calogero-Moser-Sutherland model. 
To clarify its relation to the known integrable models we have implemented two different 
types of reduction: due to discrete symmetry and due to continuous symmetry. In both 
cases we derived IIA„ Euler-Calogero-Moser-Sutherland model. Concerning the Singular 



14 



orbit strata, it was shown that in this case the corresponding dynamical system is the 
certain deformation of the Calogero-Moser-Sutherland model in a sense of description of 
particles with different masses. The masses of the particles are not arbitrary, their mass 
ratios depend on the degeneracy of the given Singular orbit stratum. As example was 
considered the case of GL{3, M) group, restricted to the Singular orbit stratum with two 
coinciding eigenvalues. In this case the reduced system coincides with the 11^42 Calogero- 
Moser-Sutherland model with particle mass ratio mi = 2m2. 
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